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Abstract. 

We have studied the effects of nonlinear scale evolution of the parton distribution 
t-H ' functions to charm production in pp collisions at center-of-mass energies of 5.5, 8.8 

and 14 TeV. We find that the differential charm cross section can be enhanced up to a 
factor of 4-5 at low pt- The enhancement is quite sensitive to the charm quark mass 
and the renormalization/factorization scales. 

m : 
o . 

Global fits of parton distribution functions (PDFs) have been obtained by several 
groups, such as CTEQ [HI2] an d MRST [HHUIS]. These fits, based on the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) jO] scale evolution, fit the deep inelastic 
scattering (DIS) HERA data [7j on the proton structure function F2(x,Q 2 ) at large 
interaction scales Q 2 > 10 GeV 2 and momentum fractions 0.005 very well. However, 
at small scales, Q 2 < 4 GeV 2 , and at small momentum fractions, x < 0.005, these fits are 
usually worse. In addition, the next-to-leading order (NLO) gluon distributions become 
very small or even negative in the small-x, small-Q 2 region. In this region, the gluon 
recombination terms, giving rise to nonlinear corrections to the evolution equations, 
become important. In a previous work |9 3 , the first of these nonlinear terms, calculated 
by Gribov, Levin and Ryskin (GLR) [TU] and Mueller and Qiu (MQ) jTT], were included 
in the leading order (LO) DGLAP evolution equations of gluons and sea quarks. It was 
shown that the HERA DIS F 2 (x, Q 2 ) data can be reproduced well with a new PDF set:):, 
EHKQS, obtained in Ref. 0, employing LO-DGLAP+GLRMQ evolution. 

Introducing the GLRMQ terms slows the scale evolution. At Q 2 < 10 GeV 2 
and x<0.01, where the nonlinearities are important, this gives rise to larger gluon 
distributions than in the pure DGLAP case. This is shown on the left-hand side of Fig.^ 
where we plot the scale evolution of the EHKQS and CTEQ61L [2 PDF sets for several 
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Figure 1. Left: Scale evolution of the EHKQS and CTEQ61L gluon 
distribution for various fixed values of x. The enhancement caused by the 
nonlinear terms vanishes during the evolution. Right: The gluon distribution 
function of x at Q 2 = 1.4 GeV 2 . 



fixed values of x. The enhancement, large at small scales, vanishes as the nonlinear 
terms become negligible at larger scales. Finally, at Q 2 > 10 GeV 2 , the evolution is 
clearly dominated by the DGLAP terms. On the right-hand side of Fig. Qwe show the 
gluon distributions as a function of x for Q 2 = 1.4 GeV 2 , the initial scale of the EHKQS 
set. As seen in Fig. ^ the HERA data relevant for nonlinear scale evolution suggests a 
factor of ~ 3 gluon enhancement at x ~ 10~ 4 and Q 2 = 1.4 GeV 2 . 

However, since the HERA F 2 data alone cannot distinguish between the linear 
DGLAP and nonlinear DGLAP+GLRMQ evolution, additional independent probes are 
needed. Here, we discuss how charm quark production in pp collisions could probe the 
gluon enhancement Charm production is an ideal choice since the charm mass 

is low and its production is dominated by gluons. Assuming factorization, inclusive 
differential charm cross sections at high energies can be expressed as 

dcr pp ^ C cx(Q 2 ,Vs) =J2 hi x uQ 2 ) ® fj(%2,Q 2 ) ®da ij ^ c7 : {k} (Q 2 ,x 1 ,X2) (1) 

i,j,k=q,q,g 

where aij^ C c{k}(Q 2 , ^i, x 2 ) are the perturbatively calculable partonic cross sections for 
charm production at scales Q 2 ~ ^> Aq CD , x\ and x 2 are the momentum fractions 
of the partons involved in the hard scattering and fi(x,Q 2 ) are the free proton PDFs. 
We assume that the renormalization scale and factorization scale are equal. Only the 
gg and qq channels are allowed at LO, which we consider here. 
We calculate the ratio of differential cross sections, 

d 3 a{EHKQs)/{dp T dydy 2 ) f . 

R(Pt, V, V2) = wt; , , r , 2 

a 6 a (CTEQ61L )/ (dp T dy ay 2 ) 
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Figure 2. The ratio of differential charn cross section, R(pT,y,y2), for fixed y 
and ?/2 as a function of charm quark px at = 14 TeV (top) 8.8 TeV (middle) 
and 5.5 TeV (bottom) in pp collision. The rapidities are y = yi = (solid), 
y = 2, i/2 = (dashed) and y = 2/2 = 2 (dot-dashed). 

where j>t is the charm quark transverse momentum and y and 2/2 are the rapidities of 
the charm quark and the antiquark. The results for the enhancement are plotted in 
Fig. El as a function of pt for fixed y and 2/2- The center-of-mass energy yfs is varied 
from 14 TeV (top) to 8.8 TeV (middle) and 5.5 TeV (bottom). Obviously, the largest 
enhancement is obtained at the largest ^/s where the x values are smallest. Different 
charm masses and scales are used, m c = 1.2 GeV and Q 2 = 4m\ on the left-hand side 
and m c = 1.3 GeV, Q 2 = m\ on the right-hand side Both of these scale choices 
lie within the applicability region of the PDFs. The rapidities are y = 2/2 = (solid), 
y = 2, 2/2 = (dashed) and y = 2/2 = 2 (dot-dashed). For the highest energy, 14 TeV, 
the maximum a/s at the LHC, the enhancement at pt ~ is a factor of ~ 5 for m c = 1.3 
GeV and ~ 1.5 for m c = 1.2 GeV. We repeated the calculations for larger masses, up to 
m c = 1.8 GeV for both Q 2 = m\ and Q 2 = Am 2 -,. We found smaller enhancements, ~ 2 
and ~ 1.25 at pr ~ 0, respectively. Clearly, the charm enhancement can be substantial, 
but it is very sensitive to the choice of mass and scale. It also vanishes rapidly with p?. 

Integrating over the rapidities does not change the result much. The maximum 
enhancement for the rapidity-integrated ratios at small pt at = 14 TeV is still ~ 4.5 
for m c = 1.2 GeV, Q 2 = 4m 2 , and ~ 1.3 for m c = 1.3 GeV, Q 2 = m\ [9]. The ratio 
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of the cross sections integrated over p T and y 2 is rather flat as a function of y, but is 
reduced to factors of 1.2 and 1.8 for the two low mass cases studied. 

Since the DGLAP gluon distributions are already well constrained by the HERA 
data, they cannot absorb additional large effects. Therefore we can conclude that, if this 
small-pr enhancement in the charm cross section relative to the DGLAP-based result 
is observed in the future experiments e.g. at the LHC, it is a signal of nonlinear effects 
on the PDF evolution. 

As discussed in Ref. the EHKQS PDFs were obtained in the region where the 
nonlinear terms do not dominate the scale evolution. However, at the smallest x and Q 2 
of the HERA data [Zj one is already close to the gluon saturation region, where all orders 
of nonlinearities become important [T^j. Therefore these PDFs, and consequently also 
the charm enhancement, should be considered as upper limits. Furthermore, a full NLO 
DGLAP+GLRMQ analysis would be needed for computing the charm hadroproduction 
enhancement consistently to NLO. 

In a follow-up work [H], we show that more than half of the charm enhancement 
survives hadronization to D-mesons. We also show that in the most optimistic case 
the enhancement can be observed in the D° pt spectrum in the ALICE detector at the 
LHC. 
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